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Abstract 

Massive multi-input multi-output (MIMO) can support high spectral efficiency (SE) with simple 
linear transceivers, and is expected to provide high energy efficiency (EE). In this paper, we analyze the 
EE of downlink multi-cell massive MIMO systems under spatially correlated channel model, where both 
transmit and circuit power consumptions, training overhead, channel estimation and pilot contamination 
(PC) are taken into account. We obtain the maximal EE for the systems with maximum-ratio transmission 
(MRT) and zero-forcing beamforming (ZEBE) for given number of antennas and users by optimizing 
the transmit power. The closed-form expressions of approximated optimal transmit power and maximal 
EE, and their scaling laws with the number of antennas M are derived for the systems with MRT and 
ZEBE. Our analysis shows that the maximal EE decreases with M for both systems with and without 
PC, but with different descending speeds. Eor the system without PC, the optimal transmit power should 
be configured to increases with M, while for the system with PC, the optimal transmit power should be 
configured as a constant independent from M. The analytical results are validated by simulations under 
a more realistic three-dimensional channel model. 
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I. Introduction 

Massive multi-input multi-output (MIMO) has become a promising technique for fifth-generation (5G) 
cellular networks. Simply by significantly increasing the number of antennas at the base station (BS), 
M, massive MIMO can support very high spectral efficiency (SE) with linear transceivers |[T]|. When 
M —^ oo, multi-user interference (MUI) and inter-cell interference (ICI) are randomized and averaged 
out without relying on sophisticated precoding and coordination, and the only limiting factor on improving 
SE is pilot contamination (PC) OJ-lSl. 

Energy efficiency (EE) is one of critical design goal for 5G networks 0|. To provide the same target data 
rate for a user or target sum rate for a system, MIMO systems can achieve higher EE than single antenna 
systems when only taking into account transmit power ||5l. However, MIMO technique improves the SE at 
the cost of extra radio frequency (RE) links and complicated signal processing, both of which consume 
considerable circuit power ||6l. By optimizing the configuration of spatial and frequency resources to 
maximize the EE of downlink multicarrier MIMO systems, it was shown in Q that the EE increases 
with the SE only when the frequency resource is still available, and the maximal EE of traditional MIMO 
system reduces with the increase of the number of antennas. 

Massive MIMO systems are expected to be energy efficient, since the transmit power can be signifi¬ 
cantly reduced owing to the large array gain and multi-user multiplexing gain IH. Moreover, by installing 
a large number of antennas and with the low transmit power per-antenna, inexpensive components can 
be used to build the system ||3l. In the past several years, valuable research results have been obtained 
for analyzing the EE of massive MIMO systems. Eor the uplink, it was shown in lH that massive MIMO 
has a large potential in improving EE by studying the power scaling law, where the EE is defined as 
fhe sum rafe divided by the transmit power. The results show that to achieve the same sum rate as a 
single antenna system, when M —)• oo, the transmit power of massive MIMO system reduces with the 
law of 1/M if the BS has perfect channel information and with 1/y/M if the BS has estimated channels 
with PC. When the circuit power consumption is considered, however, a massive MIMO system is more 
energy efficient than traditional MIMO systems only when the average channel gain is small or the circuit 
power consumption is low iQ. Eor the downlink, it was shown in ifTOl that the optimal transmit power of 
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BS needs to increase with M in order to maximize the EE of a massive MIMO systems with zero-forcing 
heamforming (ZEBE), with given numbers of antennas and users. The SE-EE relationship of massive 
MIMO systems was investigated in ifTTI . which showed that the EE will increase exponentially with a 
linear SE loss hy reducing the circuit power. This result can also he interpreted as: massive MIMO is not 
energy efficient to achieve very high SE, because a linear increase of SE leads to an exponential decrease 
of EE. While existing results provide useful insight on understanding the potential EE gain of massive 
MIMO, and partially answer the question of whether massive MIMO is energy efficient, a quantitative 
characterization of how the EE scales with M remains open. Moreover, the prior results are all obtained 
under the following assumptions: (1) single cell system without ICI, (2) perfect channel information or 
channel estimation without PC, and (3) independent and identically distributed (i.i.d.) channels. 

In this paper, we investigate the EE of downlink multi-cell massive MIMO systems and strive to 
answer the following questions. (1) What is the maximal EE of the massive MIMO systems with the 
given number of antennas? (2) How should the massive MIMO systems be configured in practice to 
achieve the maximal EE? (3) What is the impact of PC, channel estimation errors, multi-cell setting, 
channel correlation, and specific beamforming? 

In real-world cellular networks, the traffic load in a cell may vary with time (e.g., day and night) and 
location (e.g., urban and suburban). This implies that in general the number of users is not a parameter that 
can be configured. Therefore, we optimize the transmit power of the BS to maximize the average EE of 
the system, given maximal ratio transmission (MRT) or ZEBE together with equal power allocation among 
multiple users, which are widely accepted precoders for massive MIMO ||2][8l. To answer the proposed 
questions, we derive and analyze the closed-form expressions of the optimal transmit power and maximal 
average EE and their scaling laws with M. Since PC is the limiting factor for massive MIMO systems, 
which however can be mitigated by some recently-proposed pilot decontamination methods such as |[T2l . 
we consider both the systems with and without PC. The analytical results are validated by simulations 
with a more realistic channel model, where different-level power consumption parameters, ranging from 
those in existing systems to those that might be possible in the future, are taken into account. 

The rest of the paper is organized as follows. In section II, we provide the system and power 
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consumption models. In section III, we derive the closed-form expressions of the average data rates, 
respectively for the massive MIMO systems with MRT and ZFBF. In section IV, we find the maximal 
EE hy optimizing the transmit power and derive the scaling laws of the EE and optimal transmit power 
with respective to the number of antennas. In section V, we validate the analytical results hy simulations. 
Einally, we conclude the paper. 


II. System Model and Power Consumption Model 

A. System Model 

Consider a downlink massive MIMO system consisting of L non-coordinated cells, where each BS 
equipped with M transmit antennas serves K single-antenna users either with MRT or with ZEBE. The 
transmit power of the BS is P, which is assumed to he equally allocated to multiple users. We assume 
block fading channel, where the channels are constant in a time-frequency coherence block with T channel 
uses and independent among blocks. The channel from the l-th BS (denoted by BS;) to the fc-th user in 
the j-th cell (denoted by UEj^) is modeled as |[T3l 

hofc = Rhtjfc G (1) 


where a is the average channel gain including path loss and shadowing, xij = 1 if ( = j and xij = X 
if I ^ j with X G (0,1] reflecting the difference of channel gains from the serving BS and interfering 
BSs, R S is composed of N columns of a unitary M x M matrix with N < M, which models 

the correlated channels with N angles of arrival, p = hijk ~ CAA(0,lAr) is the i.i.d. channel vector. 


CAA(0,l7v) denotes zero-mean complex Gaussian distribution with covariance matrix Ijv, and 


Uy 


A^-dimensional identity matrix. The channel correlation matrix is = xzjopRR^ G C 


''MxM 


is the 


*In order to obtain simple form expressions that can reflect the inherent nature of the system, we consider that all users have 
the same correlation matrix as in Qs). Note that the channels among users are independent, but their instantaneous channel 
vectors are less orthogonal due to the same correlation matrix. This leads to the statistically largest MUI and ICI and therefore 
reflects the worst-case performance of massive MIMO systems. Later on, we will use simulation under more realistic spatial 
correlated channel model to valid the analytical results derived from the model in Q- 
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We consider a time-division duplex system, where the channels can he estimated with the pilot 
sequences sent hy users. During the uplink training phase, K users in the same cell transmit orthogonal 
pilot sequences each with Ttr channel uses and transmit power Ptr- The pilots received at BSj will 
he contaminated if the same set of pilot sequences used in the y-th cell are reused in other cells and 
the traditional training signal and channel estimation methods are applied. Nonetheless, if some pilot 
decontamination methods such as that in ifT^ can he applied, the PC will he largely mitigated without 
increasing the training overhead. Therefore, we consider two extreme cases for the systems with and 
without PC to simplify the analysis. 

Denote Sj as the index set of the cells using the same set of pilot sequences as the y-th cell, and 
Lpj = |5j| as the number of these cells, where UE;^ transmits the same pilot sequence as UEj^. for 
I ^ Sj. Eor simplicity, we assume that Lpj = Lp, Mj. Then: (l)/or the system with PC, Sj = {1, 2,... , L} 
and Lp = L, where the same set of orthogonal pilot sequences are reused in all the L cells; {2) for the 
system without PC, we set Sj = {j} and Lp = 1. 

The minimum mean-square error (MMSE) channel estimate can he obtained at BSj as 


\vjjk — R-jj-Q I hjjfc -f ^ h.jik + p I ’ 


( 2 ) 


where Q = (^{l+x{Lp-l))apfiK^ + Mm) \ ~ CAf(0, is the noise at the BS, and 

I G Sj\j denotes the index set of all interfering cells using the same set of orthogonal pilot sequences 
as the j-th cell. We define 7 tr = to reflect the signal to noise ratio (SNR) of the uplink training. 

^tr 

It is not hard to find fhaf hjjk ~ CAf{0,^jj), where = RjjQRjj = uapRR^, and v = 
J_ ■ The paramefer v reflecfs fhe channel esfimafion accuracy. Specifically, under more 
correlated channel (i.e., wifh larger p), wifh higher uplink fraining SNR, wifhouf PC (i.e., Lp = 1), 
wifh larger average channel gain a, and wifh weaker ICI (i.e., wifh smaller x), the value of v is larger, 
indicafing a more accurate channel esfimafe. Considering fhe orfhogonalify befween fhe MMSE channel 
estimate and estimation errors, we have hjjk = hjjk + hjjfe, where hjjk ~ CAf(0,Rjj — is fhe 
channel estimation error |[T4l . 
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The downlink signal received at UEj^ is given by 

L K 

I t' ■ 

Vjk = 


EE ^Ijk''/^Imy^lmXlm T 'k^jki (3) 

i=l m=l 

where is the beamforming vector of BS/ to UE^m, \m = ya , xim is the transmit signal 

with E{\xim\‘^} = 1, Ujk rv-' CAf{0, cr^) is the noise at the user, || • || denotes the Euclidean norm, and E{-} 


denotes the expectation. Eor MRT, Eor ZEBE, is the m- th column of 

and H;; = • • • > ^hk] is the estimated channel matrix at BS; for all the K users in the /-th cell. 

Every user experiences MUI (when using ZEBE, it comes from the channel estimation errors), coherent 
ICI caused by PC, and non-coherent ICI generated by the signals from the interfering BSs. Erom (l3]l, 
the signal to interference plus noise ratio (SINR) of UEj^ can be expressed as 

^jkWiifjk^jkf 


Ijk — 


^ XikWh^jk^ikf + Yh 


kmWWjk^l 


|2 I iVo-s > 


(4) 




{!SSj ,m} 

'-- 


coherent ICI caused by PC MUI and non-coherent ICI 

where {/ £Sj,m^k}L){l ^Sj,m} denotes the set of users transmitting different pilot sequences with 
\JEjk. 

Then, the average sum data rate of the j-th cell with transmit power P at BSj can be obtained as 

K K 

R,{P) = i? Eh {log2 (1 + 7ifc)} = 5] Rjk{P), (5) 

k=l k=l 

where B is the system bandwidth, and Rjk{P) is the average data rate of UEj^. 


B. Power Consumption Model 

The power consumed for downlink transmission by a BS and the circuit powers for operating the BS 

in transmitting and receiving phases can be modeled as Ill51ll6r 

(1 - ^ + (1 - ^^)PbB 2 + + M(PaF + (1 - J 

^BS = - — --rr:;-^-, (6) 


(1 - (Tdc)( 1 - crMs)(l - <7cool) 


appropriate power consumption model is critical for evaluating the EE of a system. According to the analysis in nn . the 
circuit power consumption of a massive MIMO BS can be modeled in the same way as a traditional BS, but the corresponding 
parameters will differ. Note that in the feeder loss is removed considering that massive MIMO systems will employ active 
antennas where the RE module is integrated into antenna, which is different from a traditional BS with passive antennas El. 
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where rjpA is the power amplifier efficiency, o■DC^ <7 ms and cJcooi are respectively the loss factors of 
direct-current to direct-current power supply, main supply and cooling ifTSl . Pbb 2 is the baseband signal 
processing power consumed for computing beamforming vectors, PcE is the signal processing power 
consumed for channel estimation, P^f is the RF power consumption, and PeBid and PbBiu are respec¬ 
tively other baseband processing power consumption in the downlink transmission phase and in the uplink 
training phase. According to ifTSl . Pbb 2 can be modeled as Pbb 2 = )Rfiops,o ^ Rfiopsfl 

is the floating-point operations per-second (flops) per-antenna for each user, r]c is the power efficiency 
of computing measured in flopsAV, and we use a binary variable 6zf to differentiate MRT {6zf = 0) 


and ZFBF {6zf = 1) in computational complexity. PcE can be modeled as Pqe = 


_ Mlog^(irTtr)J7/iopa,0 


nc 


In order to differentiate the impact of transmit and circuit power consumptions on EE, we define an 
equivalent circuit power consumed at each antenna except beamforming as Pc — ' 

and an equivalent circuit power consumed for beamforming at each antenna for each user as Pgp = 

—n- s/i -r. Then, ® can be rewriften as 


Pbs = (1 - R^)pP + MPc + M{{1-R^){K + 5 zfK^) + ^log2(Prtr))P 


sp 


transmit power 

(1-^)7?P + MPo, 


circuit power 


where q = 


T 

1 


(7) 


„Po = 


^^RF + ^^BE 


^ - (i-.oc)(i-<XMs)(^-^coo.) + (^ + )Psp is an equivalent 

circuit power consumption per-antenna, and the approximation comes from the similarity of the signal 
processing power consumed in uplink training and downlink transmission ifT^ . 


C. Downlink EE 

The downlink EE is defined as the ratio of the average downlink throughput to the total power 
consumption of the P BSs, where the throughput is the sum rate of all cells excluding the uplink training 
overhead. Erom ([5]l and (|7]l, we can express the downlink EE of the network as 


EE(P) = 


L((l-^)r?P + MPo)' 


( 8 ) 


















III. Average Data Rate Analysis 


To obtain a closed-form expression of EE(P) for optimization, we derive the asymptotic data rate for 
large system, where M and K grow infinitely while M/K is finite. According to the random-matrix 
theory |[T^ . the asymptotic rate converges in mean square to the average rate. Hence, we can use the 
asymptotic data rate as the average rate. 


A. Average Data Rate for MRT 

For the system using MRT, = ^jjk and the SINR of in ([Hi can be rewritten as 


^fk 


l£Sj\j {l&Sj,m,j^k}\j{l^Sj,m} 


(9) 


coherent ICI caused by PC MUI and non-coherent ICI 

By using the methods of asymptotic analysis in |[T3l but reserving the terms related to channel 


estimation errors, we can derive the asymptotic data rate provided by each BS with MRT as 

\ 


R^s = BKlog^{l + 


tM tM I I ’ 

' ^nP ' vaP/ 


( 10 ) 


where v = — . l i < 1 reflects the accuracy of the estimated channel, and 


SM^M + ^ + {Lp-l)p, 


7trQ: 


^ x{Lp - 1)(M - p), ^ {KL - l){pLp + 


( 11 ) 

( 12 ) 


are respectively the average receive powers of the desired signal, coherent ICI, MUI and non-coherent 
ICI, all normalized by va, Lp = Yli^Sj Mi = l + x(-^P“l) and L = Xij = l + x(-^“l)- Both the 
average signal power tiers'^ = a{{M—p)v+p) and the coherent ICI power ua/p = avx{Lp — l){M—p) 
increase with v, which means that they become higher when the channel estimate becomes accurate, and 
the non-coherent ICI power val^p = a{KL — l)p does not depend on the channel estimation accuracy. 


B. Average Data Rate for ZFBF 

For the system using ZFBF, the beamforming matrix of BSj is Wj = Hjj and the 

beamforming vector for can be expressed as = YS^Yijjk, where 11? = 
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ll20l . From (O, the instantaneous SINR of UE^fc with ZFBF can he obtained as 

Z ^_ Ajt||hgtn|h„t|p _ 

E A,t||hff nfhmIP + E A„„||h"nfh„„||2 + i^' 

l&Sj\j ^ {l&Sj,m^k}U{l(^Sj,m} ^ 

The asymptotic data rate provided hy each BS with ZFBF can he derived as (see Appendix 0 

•^BS =BKlog2 + I Xo-2 ) ’ (^^) 

where the average receive powers of the desired signal, coherent ICl and the sum power of MUl and 
non-coherent ICl normalized hy va are respectively 

5^^M + :^ + (Zp-l-iT)p, (15) 

I^^x{Lp-l)iM-2pK), I^p^(KL-l){pLp + ::^)-p{K-l). (16) 

Similar with MRT, both the signal and coherent ICl powers increase with v. However, different from 
MRT, the sum of the noncoherent ICl and MUl power val^p = a{KL — l)p — vap{K — 1) decreases 
with V, since accurate channel estimate induces less MUl. 

Remark 1 (Maximum data rate for given P and infinite M): By comparing the asymptotic rates of 
MRT and ZFBF, we can find that for arbitrary system configuration, S^>S'^, Ip >/p and I^p> I^p- 
Consequently, limM^oo <5^ = limM^cxD <5^ = M, but I^p > I^p still holds that is independent from M. 
For the system without PC (i.e., Lp = 1), we have Ip = /p = 0, and hence the data rate achieved by 
ZFBF is higher than MRT, where the rate gap approaches to a constant, which is limM^oo(^BS “ ~ 

BKlog 2 ( 1+ ) ■ For the system with PC (i.e., Lp > 1), we can obtain limM^oo(-fp= 

V Cp + i;^/ 

limAf^oo(f^P+f^^p) = x(-^P“l)-^^ which indicates the same data rate for MRT and ZFBF by recalling 
that limM^oo <S'^ = limM^oo S'^- 

Remark! (Maximum data rate for given M and infinite P): From (1101) and (1141) . we have limp^oo R^s ~ 
BK log 2 ^1 + 7 ^^ and limp^oo^is ~ BR^og 2 ^1 + for MRT and ZFBF, respectively, where 
jM _ jM and = Ip +l'^p- For the single-cell systems with perfect channel estimation (i.e., 

7tr —)• oo) considered in IHOllllH . we can obtain that = {K — l)p due to MUl for MRT and 1^ = 0 
for ZFBF, which indicates that with MRT the system is interference-limited for large P and achieves a 
finite maximum data rate, but with ZFBF the rate of the system increases with P without limit. However, 
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for the multi-cell scenario with channel estimation errors, the results are different. First, we have > 0 
and /^ > 0 in multi-cell systems, because L > 1 and Z > 1. Second, imperfect channel estimation will 
lead to > 0 and /^ > 0 even for single-cell systems. Therefore, once ICI or channel estimation 
errors is considered, both systems with MRT and ZFBF become interference-limited, whose maximum 
achievable rates are finite even when P —oo. 


IV. Maximal Energy Efficiency Analysis 


In this section, we optimize the transmit power to maximize EE, respectively for the massive MIMO 
systems with MRT and ZEBE. We then analyze how to configure fhe fransmif power fo maximize fhe 


EE for given numbers of anfennas and users and fhe resulting maximal EE. To show fhe pofenfial EE 
gain of massive MIMO, we do nol consider fhe fransmif power consfrainf. Based on ([8]l, (ITOl) and (ITdl) . 


fhe maximal EE can be found from fhe following problem 


max EE(P) 


(1_ A^)BiFlog2 + 

(1 - ^)i]P + MPo 


(17) 


s.t. P ^ 0, 

where P = in (E) and / = Ip+I^p in ^ for MRT, P = in O and in (O for 

ZEBE, and G = for bofh. 

If is nol hard lo show fhaf fhe numerator of EE(P) in (ITtI) is concave, fhe denominalor is convex, and 
bofh are differential. Hence, fhe EE is a pseudoconcave function wilh respecf to P. This suggesfs fhaf 
fhe problem has a globally opfimal solution, which can be found from fhe Karush-Kuhn-Tucker (KKT) 
condition. 

If fhe circuif power consumpfion MPq = 0, if is nol hard lo oblain from fhe KKT condition fhaf 
fhe opfimal fransmif power P* = 0, Ihen fhe corresponding dala rale is zero. The maximal EE can be 
oblained by E’Hospilal’s rule as 


EE* 


BK 


lim 

(1 - ^^)ri P^o 


log2 (l + 


SP 

IP+G 


PAriog 2 e S 


(18) 


P 
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If the circuit power consumption MPq > 0, which is true in practice, after some regular manipulations 
we can obtain the KKT condition as 


P* + 


MPn 


SG 


(1 _ / ((5 + I)P* + G){IP* + G) 


- In 1 + 


SP* 


IP* + G 


= 0 . 


(19) 


A. Optimal Transmit Power Analysis 

By expressing ^ function as P* from ([T^ . it is easy to show that its first-order derivative 

over P* is positive, which gives rise to the following proposition. 

Proposition 1: P* increases monotonically with 

That is to say, in order to maximize the EE of the massive MIMO system the transmit power should 
increase with the number of antennas M, equivalent circuit power consumption per-antenna Pq, training 
overhead KTt^, and equivalent power amplifier efficiency 

To understand the behavior of massive MIMO systems with given number of antennas, we derive the 
scaling law of the optimal transmit power and maximal EE with respect to M. To this end, we need to 
find fhe closed-form expression of P* from fhe franscendenfal equation in ([T^ . which is very difficulf if 
not possible. To tackle the difficulty, we introduce some approximations, which are accurate for massive 
MIMO systems. With the derivations in Appendix |B] the approximate optimal transmit power can be 
obtained as 




+ Lp) 


1 

S+1 




( 20 ) 


which is accurate for large value of M. 

Since massive MIMO can support high SE, a natural concern is whether the operating point of P* 
that achieves the maximal EE will cause a SE loss. To address this concern, in the sequel we analyze the 
optimal data rate achieved by P* and the maximal achievable rate of ZEBE or MRT obtained by setting 
P —oo. By substituting (l20l ) into (fTOl) and (fT4b . the optimal average data rate per BS with MRT and 
ZEBE can be obtained as a unified expression, which is 


~ BK \og2 1 + 


s 


I + B- 

1 ^ p. 


= BK log2 1 -h 


5 


1 + 


Pova 


■a+7) 


ln(l+|) 
M 
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and the maximal achievable rate of ZFBF and MRT can be obtained as, 


-RBS.max = lim BK log 2 1 + 


5 


P—)-oo 




= BK log 2 1 + 


5 


By substituting I = I^p + Ip and S = for MRT and I = I'^p + Ip and S = S'^ for ZFBF, the gap 
between the optimal data rate and maximal achievable rate per BS can be derived as 


= -RBS.max - RbS 

^ BK log2 ^ 
BK log 2 f 1 + 


(1 - / 1 i\ln(l + f) 


Pova 


5 + 1 


(1-lnM-ln/„ 


PqVO. 


InpM 


M 


, without PC, 


+ (1 + 


x(Fp-i)+^ J 


1 ln(l 


+ 


x(Fp-i)+^' 


( 21 ) 


, with PC, 


( 22 ) 


where InP = I^p for MRT and I^p = I^p for ZFBF, the first approximation comes from log 2 (l + 7 ) ~ 
log 2 ( 7 ) and is accurate for large SINR 7 , the second approximation is from S' « M for large value of 
M, and both are accurate for massive MIMO systems. 

From (l22b . it can be observed that /S.p decreases with the increase of Pq and M. When M —)• 00 , 
—)• 0. This give rises to the following proposition. 

Proposition 2: For the massive MIMO system with given numbers of antennas and users, the optimal 
average rate achieved by P* is close to the maximal average rate achieved by P —)• 00 . 

This implies that supporting the maximal EE of a massive MIMO system with given values of M and 
K by configuring fransmif power will cause a liftle loss of fhe maximal sum rale achieved by ZEBE or 


MRT. 


Remark 3: By selling L = 1 and 7 tr —)• 00 in ZEBE, it is easy to show from (fT^ that P = Ip+I^p = 0 
and from (1211) that Ap — 00 . This means that for single cell massive MIMO system with perfect channel 
and ZEBE, the optimal rate is far from the maximal achievable rate (that is infinite). By contrast, for 
multi-cell massive MIMO system with channel estimation errors but without PC, either MUI or non¬ 
coherent ICI will cause Ap —)• 0 for large M. When PC is in presence, Ap —)• 0, but with a faster 
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decreasing speed with M compared with the system without PC. In other words, in multi-cell massive 
MIMO systems, « .^BS.max- 

To quantify how the optimal transmit power P* should increase with M and compare with existing 
results miTOl, in the sequel we show its scaling law. By substituting I = I^p + /p and S = into 
(l20l) for MRT and I = I^p + Ip and 5 = 5^ into (l20l) for ZFBF, the power scaling law can he obtained 
after some regular manipulations. 

Proposition 3: For both MRT and ZFBF, the optimal transmit power scales with M as follows 


~ without pc, 

^!_ (23) 

~ 0(1), with PC, 

where InP = I^p for MRT and InP = I^p for ZFBF. 

Remark 4: For the single-cell massive MIMO system with ZFBF and perfect channel information at 
the BS, it was shown in ifTOl that the optimal transmit power grows proportional to By setting 

L = 1 in the expression in (|2^ without PC, the scaling law becomes This indicates that for the 

single-cell massive MIMO system with ZFBF, the imperfect CSI will change the scaling law. Moreover, 
Proposition 3 indicates that the optimal transmit power increases with instead of once the 

system suffers from the interference other than PC, and approaches to a constant independent of M 
once the system suffers from PC. Noting that such a power scaling law is very different from that in 
m because we consider the overall power consumption at the BS but only transmit power is taken into 
account in m. 


lim P* 

M^OO 


PoKu'^(—^+\ 


( 1 - 


' P~ltT 


)r)al„ 


PoKa^( 


P7tr' 

(1 - ^)ria 


-+i 


B. Maximal EE for MRT 

By substituting = I^p + Ip and (l20l ) into the objective function in ([T7] ). we can derive the 

approximate maximal EE of the system with MRT as 


EE 


M* 


(1 - ^)BK 


log 2 1 -h 




x{Lp-l)(M-p)+{KL-l){pLp + ^)+ 




M + c/i(M) 


where c = 


(1 - 


Po 


1 \ ]fp\ 

—— is a constant, and /i(M) = 


M Af 

_ / 7^ gM_|_jM 


ln(l+f^ 


( 24 ) 
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In what follows, we analyze the scaling law of maximal EE with M for the massive MIMO system 
without and with PC, respectively. 


1} Without PC: Eor the system without PC, Lp = 1, and Ip = 0, then limM^-oo fi{M) = J jw 


M 


and the maximal EE scales with M as 

KT. 


lim EE 

M^OO 


M* 


(1 _ P^)BK log 2 (M) - log 2 (/np) (1 - ^P)BK log 2 (M) 


M 


Po M Po 

which indicates that the maximal EE decreases with M asymptotically. 

2) With PC: Eor the system with PC, Ip > 0, the maximal EE scales with M as 

log2(l + - log2(l + x(L^-I)m ) 


, InM 


(25) 


lim EE’^* 
M—^OO 


Po 

KB, 


M 


(l-^)^^log2(l + ^^^) 1 

' M’ 


Pn 


(26) 


which indicates that the maximal EE also decreases with M asymptotically hut with a faster speed than 
the system without PC. 

C. Maximal EE for ZEBE 

By substituting 5^, = /^p + Ip and (l20b into the objective function in (fTTl) . the approximate 

maximal EE of the system with ZEBE can be derived as 


z* 


EE 


where / 2 (M) = 



^)E 

Po 

M 

M 



ln( 

psy 


M+::f^+{Lp-l-K)p 

'tr “ 

- 1 - 


x{Lp-l){M-2pK)+{KL-l){pLp+:^)-p{K-l)+ 


/ 2 (M) 


M + c/2(M) 


(27) 


With (|27]) and using similar method as analyzing MRT, we can obtain the scaling law of the maximal 
EE for the systems with ZEBE. Eor conciseness, the derivation is omitted, where the results show that 
the systems with MRT and ZEBE can achieve exactly the same scaling law. 


D. Summary of the Results 

The scaling laws of the optimal transmit power and maximal EE of the massive MIMO systems with 
respect to M are summarized in Table H 
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TABLE I 

Scaling Law of Optimal Transmit Power and Maximal EE with Number oe Antennas for MRT and ZEBE 



The results can be explained as follows. When the number of users K is given, for the system without 
PC, the optimal sum rate increases with log 2 M owing to the increased array gain, while for the system 
with PC, the optimal rate approaches to a constant independent of M owing to the coherent ICI. On the 
other hand, the overall circuit power consumption increases with M. As a consequent, the maximal EE 
of massive MIMO system with MRT or ZEBE finally decreases with M. Due to the training overhead 
that depending on K, the maximal EE does not simply increase with the multiplexing gain of K. In fact, 
as have been analyzed in II101I18II . there exists an optimal number of users for fixed values of M and P 
and a fixed value of M/K. 

Based on fhe table and previous analysis, we can answer the proposed questions as follows. 

(1) Given the number of users, the maximal EE decreases with the increase of M, no matter if the 
massive MIMO system is with or without PC. 

(2) Given the number of users, in order to maximize the EE, the transmit power should be configured 
to increase with M according to the laws in Table H Since for massive MIMO system the maximal 
achievable rate increases with M but the maximal EE decreases with M, meanwhile .Rgg Ri .Rss.max as 
addressed in Remark 3, an energy efficient practical system should configure M as fhe minimum number 
of antennas that can support a required average rate, which depends on the traffic load in a specific time 
and locafion. 

(3) In single-cell scenario, fhe channel esfimafion errors change fhe scaling law of P* for massive 
MIMO system with ZEBE. In multi-cell scenario, PC changes the scaling laws of P* and EE* for both 
massive MIMO systems with MRT and ZEBE. Channel correlation affects the power scaling law and the 
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impact for the systems with MRT and with ZFBF differs, hut does not affect the EE scaling law. 

V. Numerical and Simulation Results 

In this section, we validate previous analytical results via simulations. 

Consider a massive MIMO system consisting of seven macro cells each with radius 250 m, where 
one central cell is surrounded hy six cells. In each cell 10 users are located randomly with a minimum 
distance of 35 m from the BS. The path loss model is set as 35.3+37.61og^Q d dB, where d (in m) is 
the distance between the user and BS. Both the uplink and downlink noise variances are -174 dBm/Hz 
II2TI . and the system bandwidth is 20 MHz. The length of a coherence block is T = TcBc = 1500 
channel uses, where = 3.8 ms is the coherence time for 60 km/h moving speed and 2 GHz carrier and 
Be = 400 kHz is the coherence bandwidth for urban macro cell 11211 . In the uplink training phase, the 
transmit power is Ptr = 200 mW, and Ttr = 1 channel use. MMSE channel estimator is used. Since in 
the literature there are no specific power consumption parameters for a massive MIMO BS, we evaluate 
the EE of massive MIMO systems with parameters in the year 2012 and the predicated values in the 
year 2020 provided by GreenTouch consortium, which are respectively = 1.42 W and Psp = 3.1 mW 
(2012) and Pc = 0.2 W and Pgp = 0.4 mW (2020) ifT^ . and rj = 2.51 obtained from the PA efficiency 
i/pa = 50% ifT^ and loss factors (Jdc = 6 %, cJcooi = 9% and cjms = 7% for macro BS ifTSl are the 
same for both years of 2012 and 2020. The performance of the central cell is evaluated. Unless otherwise 
specified, fhese simulafion setups will be used throughout the simulations. 

A. Accuracy of the Approximations 

To show the impact of the approximations on the analysis, we compare the maximal EEs respectively 
obtained from the following approaches in Eig. [H where the channel model in o with a = , 

p = 2 and x = 0.1 is used. 

a) The optimal solution with legend “Simu Opt” is obtained through exhaustive searching P* from 
an EE maximization problem without any approximations, where the average data rate per user is 
obtained from simulation by averaging over small scale channels. 
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Fig. 1. Evaluation of the approximations, = 1.41 W and = 3.1 mW (2012), with MRT. 

b) The numerical results with legend “Num Opt” are found respectively from optimization problem 
(fTTI) for the system using MRT and ZFBF by bisection searching, where the asymptotic data rates 
are used. 

c) The approximated solutions with legend “Num SubOpt” are computed with the approximated max¬ 
imal EEs in dH for MRT and (EH) for ZEBE. 

The results show that the approximations are accurate for massive MIMO systems with MRT. The 
results for the system with ZEBE are similar and are not shown to make the figure clear. 

In order to demonstrate that the analytical results are also valid for more realistic channels, in the 
rest of the subsections we simulate the performance of massive MIMO systems by averaging over large- 
scale and small-scale channels, where a three-dimensional (3D) MIMO system with uniform rectangular 
array in urban macro (UMa) scenario (referred as 3D UMa) is employed |[2H . considering that such 
a large number of antennas may not be arranged as a linear array. Note that the lognormal distributed 
shadowing with 4 dB deviation is considered in the 3D UMa channel, and the users are randomly placed 
in each cell. In the simulation, PC will be automatically generated in the massive MIMO systems unless 
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Otherwise specified. The antenna spacing at the BS is half of the wavelength for both horizontal and 

vertical directions. The main 3D MIMO channel parameters are listed in Table |II] 

TABLE II 

Main Parameters of 3D MIMO Channels 


Mean of azimuth clusters 

(7(-60°,60°) 

Mean of elevation clusters 

[/(-45°,45°) 

Mean of log delay spread (DS) (logj^g([s])) 

-6.62 

Variance of log DS 

0.32 

Mean of log azimuth spread (AS) (logio([°])) 

1.25 

Variance of log AS 

0.42 

Mean of log elevation spread (ES) (log]^o[°]) 

max[-0.5, -2.1 -I- 0.9] 

Variance of log ES 

0.49 

Number of clusters 

12 

Number of rays per cluster 

20 


B. EE Comparison Between Massive MIMO and LTE Systems 

To validate the analytical results and compare the EE of massive MIMO with traditional MIMO systems 
(say, ETE systems), we show the EE-Rate curves under various settings. The average data rate per BS 
with MRT or ZEBE and corresponding EE of the system are obtained from simulations under 3D UMa 
channel, where the transmit power P increases from 0 W without limit and P* is found by exhaustive 
searching. To show the feasible region of the rate achievable by a macro BS with a maximal transmit 
power of 40 W (the maximal transmit power for a LTE macro BS), the results with P > 40 W are 
plotted with dotted curves. 

In Eig. |2] (a), we show the impact of the number of transmit antennas M as well as the impact 
of optimizing the transmit power and optimizing the number of users, where the power consumption 
parameters provided by GreenTouch consortium for the year 2012 are considered (the results with the 
parameters for 2020 are similar and hence are not shown). 

We first fix the number of users, considering the application scenario with low traffic load where user 
scheduling is no necessary. To be comparable with at least one kind of LTE BS, say macro BS with 
M = 8, we set K = 8. Due to the MUI, non-coherent ICI and coherent ICI, the data rates of the massive 
MIMO systems are limited to finite values when P —)• oo for both MRT and ZEBE, as addressed in 
Remark 2. The maximal EEs are achieved approximately at the maximal achievable rates of ZEBE and 
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MRT, which agrees with Remark 3. Before the sharp reduction, the EEs achieved hy MRT and ZEBE 
are very close, which can he interpreted hy Remark 1 and their similar level of circuit powers (with only 
difference in Psp)- The maximal EE achieved hy ZEBE exceeds that hy MRT due to the higher optimal 
rate, hut their maximal EEs will finally become identical (not shown in the figure) again due fo fhe facf 
in Remark 1. The maximal EEs decrease wifh M, and fhe optimal fransmif powers are much less fhan 
40 W. Eor example, for massive MIMO system wifh ZEBE, when M - 256, P* - 3.6 W. The achieved 
EE is sensifive fo fhe fransmif power, especially when if is configured heyond fhe opfimal value. 

In practice, fhere may exisf large number of users in a macro cell, e.g., in a hof-spof area during busy 
hours in a day. In such a scenario, we can furlher optimize K fogefher wifh P fo maximize EE, which is 
achievable by using an opfimal user scheduler. In Eig. |2] (a) we also give fhe EE-Rafe relafionship of fhe 
system wifh fhe opfimal number of users K* (see fhe group of curves marked by “K = K*”), where K* 
is obfained by joinfly opfimizing {K, P) fo maximize fhe average EE when fhe number of antennas M 
is given. Specifically, fhe average EE is obfained by simulation under 3D UMa channel, and {K*,P*) 
is found by exhaustive searching. If can be observed fhaf K* for fhe sysfem wifh MRT is larger fhan 
fhaf wifh ZEBE, and bofh exceed K = 8. The sysfem for a given M wifh K* achieves higher EE fhan 
fhe sysfem wifh iT = 8 as a resulf of higher rate, bofh for ZEBE and MRT, while opfimizing K brings 
higher EE gain for MRT fhan for ZEBE. When M = 256, fhe sysfem wifh MRT and K* (which equals 
fo 35) can achieve bofh higher maximal EE and higher achievable rate fhan fhe sysfem wifh ZEBE and 
K* (which equals fo 15). The maximal EE sfill decreases wifh M, which indicafes fhaf optimizing fhe 
number of users will nol change fhe scaling law of fhe maximal EE wifh M. 

TABLE III 

System Parameters oe Four Types LTE Systems ifTsll (“PL”: path loss model, Dc: cell radius) 


BS type 

System and Power Consumption Parameters 

Macro 

M-8,K-8,D^- 250m, PL - , Pmax - 40W, P, - 31.7W, P,,, - 7.8mW, rj - 3.24, crfeed - 3dB 

Micro 

M = A,K^4:,Dc = 100m, PL = 4^^, Pmax = 6.3W, Pc = 21.4W, Psp = 23.5mW, rj = 4.04, crfeed = 0 

Pico 

M = A,K = A,D, = 50m, PL = , P„ax = 0.13W, Pc = 2.6W, Psp = 2.7mW, r? = 13.72, crfeed = 0 

Femto 

M = 2, iC = 2, Pc = 30m, PL= , P„ax = 0.05W, Pc = 1.9W, Pap = 7.2mW, 7? = 21.11, crfeed = 0 
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Average Data Rate per BS (bit/s) 
(a) EE-Rate under various M 



Average Data Rate per BS (bit/s) 
(b) EE-Rate under various Pc and Psp 


Fig. 2. EE-Rate relationship under various M and Pc, Psp- The dashed and solid lines represent the results for MRT and ZEBE, 
and the maximal EEs for MRT and ZEBE are marked by lower and upper triangle, respectively. In Eig. |^(a), Pc = 1.42 W and 
Psp = 3.1 mW (for the year of 2012). In Fig. [2](b), M = 128 for the massive MIMO BS, the EEs achieved by P = 40 W 
are marked by stars for MRT and bullets for ZEBE, and the achieved EE-Rate point by four types of LTE BSs are marked 
by squares. In both Figs. [2] (a) and (b), L = 7, 3D UMa channel model is employed for the massive MIMO system and i.i.d. 
channel model is employed for LTE systems. 
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In Fig. |2] (b), we show the impact of Pc and Pgp, where the circuit power consumption parameters in 
2020 and 2012 are considered, and K is fixed as eight. As expected, the maximal EE reduces with the 
increase of the circuit power consumption. ZEBE achieves a higher maximal EE as a result of higher 
achievable rate. Along with the increase of circuit power consumption, the rate corresponding to the the 
EE-maximizing optimal transmit power moves closer to the point with infinite transmit power for both 
MRT and ZEBE, which agrees with Proposition 2. Eor comparison, the EE-Rate points achieved by four 
types of ETE systems using ZEBE transmitting at their corresponding values of P^ax are provided, where 
the parameters given in Table 17171 are for the year of 2012, and the radiated EE HI is also shown, where 
Pc = 0 and Pgp = 0 and r/ = 1. As suggested in ifTSl as well as indicated by the EE scaling laws 
(i.e., EE is inversely proportional to Pq), the EE improvement largely depends on the value of the circuit 
powers. We can see that the maximal EE of the massive MIMO system with ZEBE and predicted power 
consumption parameters in 2020 is higher than all ETE systems (e.g., 82 times higher the ETE macro 
system), and that with the parameters in 2012 is only higher than the ETE-micro and macro systems 
(e.g., 13 times higher than the ETE macro system). 

We can also obtain a target maximal EE by changing the equivalent circuit power consumption per- 
antenna Pq, from which we can see how the value of Pq should be to achieve an expected fold of 
EE improvement over traditional systems. Eor example, for the ETE macro system in 2010, the EE is 
EElte = 9.13 X 10^ bit/Joule, where the cell spectral efficiency is 6.16 bps/Hz/cell ll^ and fhe power 
consumption is 1350 W ifTSl . If a massive MIMO system wifh 256 or 512 antennas serving 10 users 
wifhouf PC is expected to achieve 1000•EELTE^ the required value of Pq should be 28 mW (if M = 256) 
or 19 mW (if M = 512) for MRT and 37 mW (if M = 256) or 23 mW (if M = 512) for ZEBE, which 
are much smaller than the predicted circuit power per-antenna 204 mW for MRT and 244 mW for ZEBE 
in 2020 by GreenTouch ifT^ . 

C. Impact of Number of Antennas on Optimal Transmit Power 

The optimal transmit power versus M are given in Eig. [3] We simulate the optimal transmit power in five 
scenarios, which can observe fhe impacts of channel correlation, channel estimation errors, non-coherent 
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ICI and coherent ICI, respectively. The circuit power parameters for the year of 2012 are considered, 
while the results of the system with the values for 2020 are similar. 

The scaling laws of the system with ZFBF are given in Fig. [3] (a). We can see that the optimal transmit 
power increases with M. For the single-cell system with perfect CSI, the optimal transmit power with 
i.i.d. and 3D UMa channels are nearly the same for ZFBF, which indicates that spatially correlated channel 
has little impact on the power scaling law of the system with ZFBF. When the channel estimation errors 


or non-coherent ICI is considered, the optimal transmit power increases with M in the law of 
and increases much slower than the system with perfect channel information, which is consistent with 
Remark 4. When the coherent ICI is considered, the optimal transmit power approaches to a constant for 
large M. As expected, the optimal transmit power decreases when the system suffers more interference. 
The scaling laws in different scenarios are consistent with Proposition 3. 

The scaling laws of the system with MRT are given in Fig. [3](h). The optimal transmit power increases 


with M proportional to ^ for all scenarios except with PC, which approaches a constant. In contrast 
to the results for ZFBF, the optimal transmit power of the system with i.i.d channel is much smaller than 
that with 3D UMa channel, which means that the spatially correlated channel has large impact on the 
performance of the system with MRT. This is because the system under the 3D UMa channel suffers 
less interference than that under the i.i.d. channel. 

Comparing Figs. [3] (a) and (h), we can find that the optimal transmit power of the system with ZFBF 
is much higher than that with MRT. 


D. Impact of Number of Antennas on Maximal EE 

The maximal EEs of the system with ZEBE and MRT versus M are given in Eig. IH where the power 
consumption parameters in 2012 are used (the results of the system with the power consumption in 2020 
are similar). It can he observed that the maximal EE decreases with M and is proportional to for 

the system without PC and for the system with PC for large M, i.e., the maximal EE of the system 
with PC decreases with M faster than that without PC. Moreover, the system without PC is more energy 
efficient than the system with PC. These results are consistent with the previous analytical results. 
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Fig. 3. Optimal transmit powers of massive MIMO systems with ZFBF and MRT versus M, the legends “perCSF’ and “estCSI” 
respectively denote the results with perfect channel information and estimated channels. 
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(a) ZFBF (b) MRT 


Fig. 4. Maximal EEs of massive MIMO systems with ZFBF and MRT versus M, 3D UMa channel model. 

























24 


VI. Conclusions 


In this paper, we investigated the potential of downlink multicell massive MIMO system in improving 
EE hy optimizing the transmit power. We derived the closed-form expressions of approximated maximal 
EEs of the massive MIMO systems with MRT and ZEBE under a spatially correlated channel model, 


and then analyzed the scaling laws of the optimal transmit power and maximal EE with the number of 


transmit antennas M for the system with and without pilot contamination. Analytical results showed that 
the maximal EE of massive MIMO systems reduces with M, and the decreasing speeds are different for 
the systems with and without pilot contamination. The optimal transmit power should he configured to 
increase with M in order to maximize the EE of the system without pilot contamination, hut he configured 
as a consfanf independenf from M for fhe sysfem wifh pilof confaminafion. Channel correlafion has large 
impacf on fhe power scaling law for massive MIMO wifh MRT, and has minor impacf on fhaf wifh ZEBE. 
Eor practical systems wifh given traffic loads, fhe mosf energy efficienf massive MIMO system should 
he configured fo use fhe minimum number of antennas fhaf can supporf fhe required average sum rate. 
Simulation resulfs validafed fhe analytical analysis, and showed fhaf fhe conclusions are also valid under 
more realisfic channel model. 

Appendix A 

Derivation of Asymptotic Data Rate with ZEBE 

In order fo derive fhe asympfolic dafa rate wifh ZEBE, we firsl derive fhe asymplofic rale of fhe 
regularized zero-forcing (RZE) beamforming. The beamforming malrix of BSj wifh RZE is 



iMxA 


(A.l) 


where (/? is fhe regularized parameter, Ilj = -|- salisfying 11 

Iljhjjfc. The receive SINR of UEj^ wifh RZE is given by 



= Ilj, and Wjfc 



^jk\\ 


(A.2) 
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According to Theorem 3 in ll20l . the asymptotic SINR of UEj^ with RZF, 7 ^, converges to the 
asymptotic SINR with ZFBF, 7 ^^^, when the regularized parameter approaches zero, i.e., 

In the sequel, we will first derive the asymptotic SINR with RZF, and then obtain the asymptotic SINR 
with ZFBF hy letting 0. 

1) Normalizing Parameter: According to the proof of Theorem 2 and Theorem 3 in EOl . the normal¬ 
izing parameter Xjk of ZFBF can he obtained as 




- 

M,K^oq 


> M^6°, 






(A.3) 


where the symbol a.s. stands for almost surely convergence, 6j{ip) = -^ti{^jjTj{ip)), Tj{ip) = 


(f y(4») = = Ti(^) (^I + f j lAv). = 

uapRR^, and <5° = 

2) Signal Power: According to the proof of Theorem 6 in ifT^ . the power of desired signal can be 
derived as 


TT fi 




(A.4) 


where the term in the denominator asymptotically approaches to 

With the expression of hjjk given in El** the numerator of (IA.4I) becomes 

IjjhjjiJJjkhjjkl 


(A.S) 


(A. 6 ) 


+ M^^fjk^jk'RjjQ + ^^fk) 


^jjk 


+ \ hjjfcFIjfcRjjQ I hjZfc + ) hjjfcQRjjlljfch 


ijjfc f ) 


where the third term on the right-hand side converges to zero almost surely by Femma 4 (iii) in ifT^ . By 
using similar way to derive (IA.51) . we can show that the first term converges to | 




M,K^oo 


By using Femma 4 (ii) in |[T3l twice and noting that ^jj = vRjj, the second term almost 


surely converges to ^iw4^tr(RjjT'' {ip)) = 7^(5 • (</?), where T - ((/?) = Tj((/?) ( -f § 


(!-•») 




vM^ j 




-33 -r 


Tj((/9), 
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and 5 ;{^) = (y.)). Define = lim^^o {ip) = (^1 + ^ and = 

PK 1 
M-pi^ ■' 


itr($jjT"°). Upon substituting ^jj = uapRR^ and T° = (^g^RR^ + I) we can derive that 


6j = vap5^. 

Based on the above results, we can obtain the signal power as 


I M hjjfc I 


-)■ 


1 M{8,{p))‘^ + {l-l)6'Up) 


M,K^oo M 




vaiM + ::^+p{Lp-l-K)) 


(A.7) 


tp—>0 M 

3) Coherent ICI Caused by PC: By using the same way to derive the signal power and considering 
that = HiiQUij = x^ih we can show that the coherent ICI caused by PC converges almost surely 
to 

2 „ „ 1 M{x6i{p))‘^ + (^ - X)x5'i{‘p) 


A 


Ik 




M,K^oo M 




Xva{Mx + ^ + p{Lp -X- XK)) 


(A.8) 


<p^o M 

4) MUI and Non-coherent ICI: According to the proof of Theorem 6 in |[T3l . the interference caused 
by I e Sj and m ^ k of RZF is 

tr(R;,-T;'M) , , \tr:{^i^T,{ip))\^';{p) 


^lm\M^ljk^l^llml 


a.s. 1 
^Xi 


M,X—>00 


M2 




+ 


M3 


{I + 6i{p))‘^6'i{p) 


(A.9) 


_ i-mp 

M3*we 


tr:i^l,T';{ip))itr{^l,Ti{p))y 


{I + 6i{p))6'i{p) 

For I = j and m y k, when y? —)• 0, the MUI of ZFBF converges almost surely to 

A- n-h ■■ P ' {f-v)ap 

M jjk^ J 


(A. 10) 


M,K—^oo,ip —>0 M 

For I £ Sj\j and m y k, when y? —)• 0, the non-coherent ICI of ZFBF caused by the users in the cells 
using the same set of pilot sequences converges almost surely to 


^lm\ M^ljk^t^ttrnl 


(1 - Xv)xoip 


(A. 11) 


M,K—¥^oo,ip—^0 M 

For I ^ Sj and arbitrary m, the non-coherent ICI of ZFBF caused by the users in the cells using 
different set of pilot sequences can be derived as 

, Xap 


^lm\ M^ljk^l^llm 


2 a.s. ^ X [P) 


M. 


,K->oo vM 5'i{p>) <p^o M 


(A. 12) 
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By substituting (IA.7I) . (IA.81) . (lA.lOl) . (lA.l 11) and (IA.121) into (IA.2I) . we can obtain the asymptotic SINR 
with ZFBF as 


^ M,K^oo 




M + 


■ytrOl 


+ p{Lp -1-K) 


x{Lp - 1)(M - 2pK) + {KL - l){pLp + ^)- p{K-l) + ^^^ ’ 


vaP 


where Lp = 1 + xi^P — L = 1 + x{L — 1) and v =- — 


^pLp + ^ ■ 

'tr 

Then, the asymptotic data rate per BS with ZFBF can be obtained as 


^BS = = BK log2 (l + x%) . 

M,K^oo ^ j ' 


(A. 13) 


Appendix B 

Approximated Optimal Transmit Power P* 
The KKT condition of P* in (fT^ can be expressed as 

SP* 


soiP- + j^;) 

{{S+ I)P*+G){IP*+G) \^IP*+G 


MPo 

H'-l'tr ' 


G 


P* _i_ G 
^ ^ S+I 

MPo G 


P* + 

MPo 


G 

S+I 


2G 



I 

p* + f 


MPo 

G 

(1-^)7 

I 

P* + f 


MPo 

T-. N 

2G 


-In 1 + 


G 


{S + I)P^ 


+ln 11 + 


G 


V IP^ 


-Mi + f 


G 

(S+I) 


P* + 


G 


■ + 


G 

I 


2(S+I) 


-In (1 + 1 

P* + ^ \ I 


G 


P* _|_ 

^ ^ S+I 

MPo MPo 


P* + 


-In (1 + ^ 


W (1-^)^ (1- 


7777 


P* j_ ^ P* _L 

^ -P 5+7 ^ -r T 


l^-lnfl + f ) =0, 


/ 


(B.l) 


where approximation (a) comes from the first order approximation of ln(l + x) « which is accurate 
for small x. When M is very large, both S and P* are large according to Proposition 1. Furthermore, if 
the system is pilot contaminated, I increases with M and is also very large. Therefore, approximation 
(a) is accurate for large value of M. Approximation (b) comes from P* + 2 {^+i) P* + (^J+i) 


P* + ^ ^ P* + Y’ since P* )$> and P* (f for large P* and S, which are true for large value 


G 


G 


of M. Approximation (c) comes from 


MPo 


G 


^ I > S+I 


because of large values of M and S in 


massive MIMO. Furthermore, if Pq is large, the value of P* will be large, then the approximations are 


more accurate. 
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The approximate solution of (O (i.e., the solution of (IB.lll l is 



ln(l + f)’ 


(B.2) 


T /I ^ 2 tr \ 1 

where approximation (d) comes from ignoring \ because ,<j « T h-, ^ j for large values of S 

and MPq due to large M, and approximation (e) comes from ignoring the second term 5 + 5 ^^ 

because the first term (proportional to is much larger than the second term (proportional to 

According to the analysis above, approximations (a)-(e) are all accurate when the value of M is large 
and become even more accurate when the value of Pq is large. 
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